Doping graphene by adsorption of polar molecules at the oxidized zigzag edges 
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We have theoretically investigated the electronic and magnetic properties of graphene whose zigzag 
edges are oxidized. The alteration of these properties by adsorption of H2O and NH3 molecules have 
been considered. It was found that the adsorbed molecules form a cluster along the oxidized zigzag 
edges of graphene due to interaction with the electro-negative oxygen. Graphene tends to donate 
a charge to the adsorbates through the oxygen atoms and the efficiency of donation depends on 
the intermolecular distance and on the location of the adsorbed molecules relative to the plane of 
graphene. It was found that by appropriate selection of the adsorbates, a controllable and gradual 
growth of p-doping in graphene with a variety of adsorbed molecules can be achieved. 

PACS numbers: 73.22.-f,73.20.Hb,73.21.La,75.25.+z 
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I. INTRODUCTION 

The edges of graphene, in particular the zigzag edges, 
play the most important role in defining the electronic 
properties of graphene^. Geometry of the zigzag edges 
localizes the electrons with maximum of the electron den- 
sity at the border carbon atoms. The localized states 
form the flat conduction and valence bands near the 
Dirac points and the size of the gap between them is 
denned by the contribution of the confinement effect and 
shape of the edgesS^. Any changes in the edge geometry, 
such as termination, reconstruction or distortion, would 
modify the electronic properties of graphene. Termina- 
tion of the edges most often is realized through the edge 
hydrogenation, which alters the sp 2 network at the zigzag 
edgesi For the pure graphene edges, there is a mixture 
of the sp and sp 2 hybridization, while for hydrogenated 
graphene the sp 2 and sp 3 hybridization dominates pro- 
viding an increase in the size of the gap. Additionally 
to hydrogenation, structural changes in the zigzag edges, 
such as reconstruction and aromaticity, also have effect 
on the electronic properties of graphene^. Moreover, the 
geometry of pure zigzag edges are found to be metastable. 
It was shown that the planar reconstruction of the zigzag 
edge to reczag type is possible because it lowers the 
edge energy^. The electronic structure of graphene with 
reczag edges differs significantly from that with zigzag 
edges, namely by the appearance of the degeneracy of 
conduction and valence bands for the momentum ranged, 
k < 27r/3, and their dispersion for k > 2ir/3 (the zigzag 
edges imposes the flat degenerate conduction and valence 
bands within the interval^, 2ir/3 <| k |< it). Coexistence 
of the zigzag and reczag edges in graphene samples has 
been confirmed experimentally 6 . 

In fabrication of graphene-based devices the control 
of the shape and quality of the edges is clearly an im- 
portant issue. With the standard approaches devel- 
oped to create graphene nanoribbons, such as litho- 
graphic patterning 7 , chemical vapor deposition method 8 
and chemical sonication^, the control of the nanorib- 
bon size and quality of the edges is poor. Most recent 
achievement in fabrication of graphene nanoribbons is 



by unzipping the carbon nanotube o 10 ' 11 . Unzipping of 
multiwalled or singlewalled carbon nanotubes by plasma 
etching provides smooth edges and a small range of the 
nanoribbon width 1 - (10-20 nm). The second unzipping 
method - the solution-based oxidative process allows to 
perform the longitudinal cut of the nanotubes and to 
obtain the nanoscale structures characterized by the pre- 
dominantly straight linear edg However, because of 
the applied solution-based oxidative method, the zigzag 
edges of graphene are oxidized and so the surface of 
graphene can be partially oxidized. Annealing or chem- 
ical reduction of samples in the N2H4 environment re- 
duces the number of oxygen-containing functionals at the 
surface and edges of graphene and, indeed, a significant 
increase in conductivity of graphene has been obtained 11 . 

From our point of view, an important improvement 
of oxidation of the zigzag edges of graphene is the pre- 
vention of reconstruction of the zigzag edges. However, 
an extensive investigation of the electronic properties of 
graphene with oxidized zigzag edges is not yet avail- 
able in the literature. Therefore, the present study is 
devoted to this particular issue and the possibility of 
doping nanoscale graphene through adsorption of wa- 
ter and gas molecules. The spin-polarized density func- 
tional theory with semilocal gradient corrected functional 
(UB3LYP/6-31G3) in the Jaguar 6.5 program^ has 
been applied in our work. The impact of the van der 
Waals interactions is not included within the DFT that 
leads to underestimation of the adsorption energy^, but 
it should not undermine the reliability of DFT for inves- 
tigation of the influence of adsorption on the electronic 
properties of graphene. 



II. UNIQUE PROPERTIES OF GRAPHENE 
WITH OXIDIZED ZIGZAG EDGES 



For graphene of finite size with pristine edges, the ef- 
fect of localization of the electron density at the zigzag 
edges is a result of specific boundary conditions for the 
wavefunctions. In particular, because all border carbon 
atoms across a single zigzag edge belong to the same 
sublattice, the wavefunctions at the zigzag edge vanish 
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(a) a- spin state 



(b) a- spin state 



(c) spin density 
charge charge 




FIG. 1: (color online) Electron and spin density distribution in a pure graphene flake and graphene with oxidized zigzag edges. 
The electron and spin densities are plotted for isovalues of ±0.03 e/A 3 . (a) The electron density distribution in the pristine 
graphene for a state of C21, symmetry (a-spin for HOMO), (b) The electron density distribution in graphene with oxidized 
zigzag edges for a state of Ci symmetry (a-spin for HOMO-9). (c) The spin density distribution in graphene with oxidized 
zigzag edges for a state of Ci symmetry. Localization of spins on the oxygen atoms obtained within the NBO analysis and the 
natural charge for the oxygen and carbon atoms within the C-0 bonds are also indicated. 



on a single sublattice, while at the armchair edges on 
both sublattices. The generated localized states at the 
zigzag edges are responsible for formation of the conduc- 
tion and valence bands near the Dirac points. Depending 
on the spin ordering of the localized states between the 
zigzag edges, graphene of finite size can be gapless or has 
a gap. If spins of the localized electrons are antiferro- 
magnetically ordered between the opposite zigzag edges 
for which the border carbon atoms belong to different 
sublattices, the sublattice symmetry is preserved. In the 
case of an absence of the confinement effect, graphene 
with preserved sublattice symmetry is expected to be 
gapless. The ferromagnetic ordering between the zigzag 
edges leads to breaking of the sublattice symmetry and 
opening of a gap. In the case of broken sublattice sym- 
metry (for nanoscale graphene the sublattice symmetry is 
broken already for C2„ symmetry when the mirror plane 
of symmetry is perpendicular to the zigzag edge s 15 ' 16 ), 
the a-spin state is localized on one zigzag edge, while the 
/3-spin state on the opposite zigzag edge. We presented 
an example of the distribution of the electron density of 
the localized a-spin state for the highest occupied molec- 
ular orbital (HOMO) in Fig.[T](a). A similar distribution 
but with localization of the electron density on the op- 
posite zigzag edge is observed for the lowest unoccupied 
molecular orbital (LUMO). The maximum of the electron 
density of these states is located on the border carbon 
atoms and, as a result, the spin density is localized there 
as well (see Fig. 2 in Refi^ 5 -). 

Oxidation of the zigzag edges makes a difference in 
the electron density distribution of the localized states 
because of alteration of the boundary conditions for the 
wavefunctions. The electron density is mostly localized 
on the sides of the polar C-0 bonds, as it is presented 
in Fig. [T] (b). For the structure shown in Fig. [1] (b), 
the electron density for the HOMO and LUMO are de- 
localized over the whole graphene structure, while the 



localized states are shifted deeper into the bands from 
the gap edges, i.e. to the HOMO-9 and HOMO-10 in 
the valence band and to the LUMO+1 and LUMO+2 in 
the conduction band. However, an increase in the size of 
the structure leads to backward shifting of these states 
closer to the band gap. Therefore, when we increased the 
width of nanoscale graphene by adding two rows of the 
carbon hexagons along the single zigzag edge, the local- 
ized states in the valence band has already been shifted 
to HOMO-5 and HOMO-6 orbitals. 

The spin density for graphene with oxidized edges is 
mostly localized on the oxygen atoms, as it is presented 
for nanoscale graphene in the Ci symmetry state in Fig. [I] 
(c) . Ordering of spins of the localized electrons along the 
zigzag edges and between opposite edges is defined by the 
applied symmetry. For the Z?2h symmetry, the spin or- 
dering is ferromagnetic along the edges and between the 
edges, while for the Civ symmetry (when the symmetry 
plane is perpendicular to the zigzag edges) and for the C\ 
symmetry, there are mixed spin ordering along the zigzag 
edges and antiferromagnetic ordering between the zigzag 
edges. The total energy of the system decreases with low- 
ering of the symmetry of graphene and the C± symmetry 
state is characterized by the lower energy. For the state 
of the C\ symmetry, the total spin at each zigzag edge is 
nonzero, particularly the total spin equals +1.0 for one 
edge and -1.0 for another edge. Because the spin density 
is mostly localized on the oxygen atoms within the C-0 
bonds, it is negligible for all carbon atoms in graphene. 
Symmetry also has influence on the size of the band gap. 
For the size of graphene presented in Fig. 1 (b), the gaps 
are AD 2h =0.18 eV, A Cjh) =1.0 eV and A Cl =1.33 eV. We 
consider that graphene obtained by nanotube unzipping 
would have highly tensile structure for which the high 
symmetry state would be unreachable, and therefore, we 
have used the Ci symmetry for all our calculations. 

We have performed the natural bond analysis 17 (NBO) 
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for nanoscale graphene with the oxidized edges and de- 
fined the natural charges for the carbon and oxygen 
atoms within the C-0 bonds, whose magnitudes are indi- 
cated in Fig. [1] (c). The C-0 is the polar bond for which 
the decentralization of the electron clouds between the 
carbon and oxygen atoms leaves the oxygen atoms neg- 
atively charged, while the carbon atoms are positively 
charged. As a result, the electro-negative oxygen atoms 
within the polar C-0 bonds are more chemically active 
than on the non-polar surface of graphene. It was found 
that for the C-0 bonds the natural charge distribution 
(see Fig.[T](c)) and the length of the bonds deviate along 
the zigzag edges (from 1.26 to 1.31 A). However, despite 
the dispersion of the bond length the variation of the 
dipole moment among the C-0 bonds is consistent with 
the deviation of the natural charge distribution, i.e. for 
the bonds characterized by a larger charge difference be- 
tween the oxygen and carbon atoms, the dipole moment 
is slightly larger as well. Therefore, interaction of the 
oxygen atoms with the adsorbed molecules would also 
slightly deviate along the edges. 



III. ADSORPTION OF H 2 AND NH 3 
MOLECULES 

The most popular and widely investigated method of 
doping graphene is by adsorption of the gas molecules 
which is claimed to be controllable by the nature and con- 
centration of the adsorbate a 18 ' 19 i 20 ' 21 i 22 ' 23 i 24 . However, 
according to the theoretical predictions, interaction be- 
tween the surface of graphene and most of the molecules 
adsorbed on graphene should be lowi* and doping of ideal 
pristine graphene (sp 2 network) should not be efficient. 
Moreover, because of weak interaction of the adsorbates 
with the ideal graphene surface the type of doping and 
its efficiency is defined by the orientation of the adsorbed 
molecules on the surfac o 15 ' 25 > 26 ' 27 . Therefore, doping ob- 
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served in experiments 
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is most probably 



the result of lattice defect in graphene^. In this respect 
we would like to emphasize two important advantages of 
oxidation of the zigzag edges of graphene: potentials to 
achieve a controllable doping (because oxidation will en- 
hance interaction of the edges with the surrounding envi- 
ronment) and stabilization of the edge geometry. The ox- 
idized zigzag edges attract polar molecules from the envi- 
ronment more strongly than the graphene surface due to 
polarity of the C-0 bonds. The charge exchange between 
the C-0 bonds and the adsorbates would contribute to 
the alteration of the electronic balance in graphene lead- 
ing to its doping. Therefore, adsorption would influence 
the graphene conductivity and offers an unique opportu- 
nity for its application in graphene-based gas sensor a 18 ' 19 . 

The H2O or NH3 molecules have been chosen for our 
investigation because of their polarity and the presence 
of more than two hydrogen atoms in their structure. Ad- 
sorption of H2 O or NH3 molecules on graphene with ox- 
idized edges has been performed for primary optimized 
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FIG. 2: (color online) Structure of nanoscale graphene with 
oxidized edges where the adsorbed molecules are attached 
to the zigzag edge by hydrogen bonds. Several adsorbed 
molecules form a cluster due to the hydrogen bonding with 
zigzag edge of the oxidized graphene and between each other, 
(a) adsorbate is the H2O molecule, (b) adsorbate is the NH3 
molecule. 



geometry of graphene. After adsorption an additional op- 
timization procedure has been applied. During this opti- 
mization the positions of the carbon atoms belonging to 
graphene are kept fixed, while the oxygen atoms within 
the C-0 bonds and the adsorbed molecules were relaxed. 
The oxygen atoms within the C-0 bonds make hydro- 
gen bonds with the adsorbed molecules due to strong 
polarity of the C-0 bonds. If adsorbed molecules have in 
their structure several hydrogen atoms and the electro- 
negative atoms possess the lone pair (oxygen, nitrogen), 
they have the ability to bind to each other through the 
hydrogen bonds, thereby forming a cluster. However, 
the hydrogen bonding with graphene leads to aligning 
of the cluster structure along the zigzag edge, as pre- 
sented in Fig. [2] (a) and Fig. [2] (b) for the H 2 and NH 3 
molecules, respectively. We have estimated the binding 
energy required to dissociate graphene from the adsorbed 
molecules. It was found that the binding energy grows 
almost linearly with increasing number of the adsorbed 
molecules, but have a slight tendency for saturation when 
the number of molecules is larger than eight. 

The highest symmetry applicable for graphene with 
adsorbates is the C\ symmetry because the charge ex- 
change with the adsorbates generates an imbalance in 
the electron distribution between the two sublattices in 
graphene and, therefore, breaking the sublattice symme- 
try. The size of the bandgap after adsorption is simi- 
lar to the value of that for pristine graphene in the C\ 
symmetry state (Ac^l.33 eV). However, if the HOMO 
or LUMO are localized on the adsorbed molecules, the 
size of the band gap is altered. Interaction of graphene 
with the adsorbed molecules induces a change in order- 
ing of the orbitals, particularly the states localized on 
the zigzag edges now belong to HOMO (HOMO-1) and 
LUMO (LUMO+1) orbitals while without adsorbates the 
localized states are shifted deeper into the conduction 
and valence bands (see discussion in Sec. Hi]) . 

Formation of hydrogen bonds between the adsorbed 
molecules and the oxygen atoms belonging to the C-0 
bonds of graphene leads to decentralization of the elec- 
tron clouds between the electro-negative oxygen from the 
C-0 bonds and the electro-positive hydrogens belong- 
ing to the adsorbed molecules. Therefore, the interac- 



4 



TABLE I: Influence of adsorption of H2O or NH3 molecules on 
spin and charge distributions in graphene with oxidized edges. 
All data are obtained within the NBO analysis. Qh 2 o and 
Qnh 3 are the charge accumulated on the adsorbed molecules, 
Sh 2 o and Snh 3 are the total spin on the cluster formed by 
the adsorbed molecules, So is the sum of spin located on the 
oxygen atoms belonging to the polar C-O bonds of graphene, 
Qtr is the total charge transfer from graphene to the adsor- 
bates. "Left" and "Right" refer to to the left and right zigzag 
edges as shown in Fig. [Tj (b). 
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Right 


Left 


Right 


Left 


Right 




1 


0.028 




0.00 




0.00 


0.00 


0.006 


2 


0.317 




-0.46 




-0.47 


1.01 


0.028 


3 


0.841 




0.48 




-1.59 


1.01 


0.084 


4 


0.856 




-1.02 




2.01 


-1.01 


0.068 


5 


0.869 




1.02 




-2.01 


1.00 


0.096 


6 


0.889 




1.00 




0.01 


-1.00 


0.073 


7 


0.882 


-0.006 


1.00 


-0.01 


0.01 


-1.00 


0.096 


8 


0.890 


0.053 


1.00 


-0.07 


0.01 


-0.93 


0.098 


9 


0.881 


0.348 


-1.00 


0.40 


0.04 


0.24 


0.132 


10 


0.888 


0.854 


1.01 


0.99 


-1.97 


0.01 


0.161 


11 


0.888 


0.864 


1.01 


1.01 


0.01 


-1.97 


0.183 


12 


0.889 


0.853 


1.00 


1.01 


0.01 


-1.95 


0.191 
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Left 
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Left 


Right 
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0.247 




-0.26 




-0.74 


1.00 


-0.018 


2 


0.906 




1.00 




-0.01 


-1.03 


0.035 
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0.922 




1.00 




0.01 


-1.00 


0.057 


4 


0.925 




1.00 




0.01 


-1.00 


0.052 


5 


0.933 




1.00 




-2.00 


1.00 


0.043 


6 


0.926 




1.00 




0.01 


-1.01 


0.048 


7 


0.925 


0.271 


-1.00 


0.28 


-0.01 


0.72 


0.054 


8 


0.921 


0.903 


-1.00 


1.00 


0.01 


-0.01 


0.127 


9 


0.921 


0.922 


-1.00 


1.00 


-0.01 


0.01 


0.109 


10 


0.924 


0.921 


1.00 


-1.00 


-0.01 


0.01 


0.100 


11 


0.924 


0.928 


-1.00 


-1.00 


1.99 


-0.02 


0.094 


12 


0.924 


0.924 


-1.00 


1.00 


-0.01 


0.01 


0.097 



tion of the adsorbates and graphene triggers a charge 
redistribution between them. Accumulation of the natu- 
ral charge on the adsorbed molecules determined within 
the NBO analysis is presented for H 2 (Qh 2 o) an d NH 3 
(Qnh 3 ) molecules in Table Q] The charge exchange be- 
tween graphene and the adsorbates has been also defined 
within the NBO procedure as Q tr — Qda — Qad, where 
Qda is the charge transfer from graphene to the adsor- 
bates and Qad is the charge transfer from the adsorbates 
to graphene. The Qda and Qad parts have been calcu- 
lated as the sum of the f2j — > f2* charge transfers between 
the donor orbital f2j and the acceptor orbital f£* defined 
within the NBO procedure as: 

Qda = J2 Q*J = E fc^/fe " e if 

i,j i,j 

where g$ is the donor orbital occupancy, ej, ej are the 
orbital energies, Fij is the off-diagonal element. The 




FIG. 3: (color online) The spin distribution in graphene with 
eleven adsorbed NH3 molecules plotted for isovalues of ±0.01 
e/A 3 . The left edge has a total spin 2.0, while for the right 
edge the total spin is zero. The adsorbed molecules maintain 
spin of -1.0 from each side of graphene, however the balance 
of the total spin between two sides of graphene with the ad- 
sorbates is maintained, i.e., -1.0 for one side and 1.0 for other 
side. 



charge occupancy transfer is considered for stabilizing or- 
bital interactions, i.e. when the second order interaction 
energy Al^j = — 2i^ 2 -/(e, — ej) has a positive sign. The 
calculated results for the charge exchange (Qtr) between 
graphene and the adsorbates (H2O or NH3 molecules) are 
also presented in Table [J The total spin of the states lo- 
calized on the oxygen atoms within the C-0 bonds (So) 
and on the adsorbed molecules (Sh 2 o or Snh 3 ) are also 
indicated in this table separately for left and right zigzag 
edges of graphene (see notations for the edges in Fig. [1] 
(b)). 

According to our studies, both graphene and the clus- 
ter formed by the adsorbed molecules are spin-polarized 
in most cases. For spin-polarized graphene, each zigzag 
edge is characterized by the non-zero total spin whose 
sign is different for the opposite edges. The sign and 
the total spin associated with the adsorbed molecules 
varies with the number of molecules. However, the bal- 
ance of the spin density between the left and right sides 
of the whole system which includes graphene with ad- 
sorbed molecules, is preserved (except for the case when 
a single water molecule is adsorbed). Thus, taking into 
account the spin localized on the zigzag edge and the ad- 
sorbed molecules (see the magnitude of Snh 3 + So and 
Sh 2 o + So in Table [TJ> the total spin +1.0 for one side of 
the system and -1.0 for other side are retained. The spin 
distribution for the oxidized graphene with eleven NH3 
molecules bound to the edges by the hydrogen bonds is 
presented in Fig. [3J Only one or two molecules within 
the cluster have the non-zero total spin and are spin- 
polarized, while the rest of the molecules are unpolarized. 
For adsorbed molecules characterized by zero spin den- 
sity the magnitude of the natural charge is almost zero 
(~ 0.05 e), but for a case of non-zero spin density the 
charge is non-zero (~ 0.5 e for each of the spin-polarized 
molecule). The dependence of spin-polarization and re- 
sulting nonzero magnetic moment on the number and 
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location of the adsorbed molecules can be useful for de- 
veloping spin- valve devices^. 

The charge transfer from pristine graphene to water 
is not expected to be efficient and was reported to be in 
the ranged, -0.02 4- 0.01 e depending on the position and 
orientation of water, while according to our calculations 
performed for the single water molecule within the NBO 
analysis, it is ~ 0.001 e (for details of interaction of water 
with pristine graphene seei^). The charge exchange be- 
tween graphene with oxidized edges and the adsorbates 
is quite efficient and increases with the number of ad- 
sorbed molecules (see Qt r in Table |T]). For adsorption 
of H 2 molecules, the charge donation from graphene 
to adsorbates increases gradually. However for a single 
H 2 molecule adsorbed at the zigzag edge, the charge 
exchange is insignificant because of the location of the 
molecule relative to the edges. A single H2O molecule 
interacts with graphene by both of its hydrogens which 
are directed to two nearest oxygen atoms at the zigzag 
edge. The resulting distance between the H2O molecule 
and the edges is large (the distance between O-H and 
O-C is ~ 2.2 A) and the H2O molecule lies parallel to 
the graphene plane. Therefore, the overlapping of or- 
bitals is weak. For two or more adsorbed molecules, one 
hydrogen of the H2O molecule makes a hydrogen bond 
with the other H2O molecule, while the second hydro- 
gen bonds to the oxygen atom within the C-0 bond of 
graphene. The cluster that formed is placed close to the 
zigzag edge (the hydrogen bond of O-H • • • O-C is ~ 1.2 
A) and water molecules are located not exactly in the 
plane of graphene but slightly below and above, hence 
providing an effective overlapping of orbitals between the 
interacting species (see Fig. [2] (a)). 

For the NH3 adsorbates, a similar behavior is obtained 
for up to three NH3 molecules. The most efficient charge 
exchange occurs when two or three molecules are placed 
along the zigzag edge, while a further increase of their 
number suppresses the charge transfer. The effect is re- 
lated to clustering of the NH3 molecules along the zigzag 
edge. Two or three NH3 molecules adsorbed along the 
zigzag edge are located slightly below and above the 
graphene plane somewhat similar to the water molecules 
in Fig. [3 (a) (the hydrogen bond of N-H ■ ■ • O-C is ~ 2.0 
A) generating an efficient charge exchange with graphene 
due to strong overlapping of their orbitals. However, ad- 
sorbates move apart from the plane of graphene when 
more than three molecules are added (as presented in 
Fig. [2] (b)) and the length of hydrogen bonds for some 
molecules increases (the hydrogen bond of N-H ■ ■ • O-C 
is ~ 2.5 A). Thus, in the cluster containing more than 
three NH3 molecules, they are located significantly out of 
the plane of graphene that lead to diminishing efficiency 
of the charge exchange of each adsorbed molecule with 
graphene. 

The charge exchange between graphene and the ad- 
sorbed molecules shifts the conduction and valence bands 
of graphene in the energy scale. In Fig. QJ we have plotted 
the orbital energies of HOMO and LUMO as a function of 
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FIG. 4: Effect of adsorption on the orbital energies of the 
HOMO and LUMO. The results are presented for both the 
a- and /3-spin states, (a) adsorbate is the H2O molecule, (b) 
adsorbate is the NH3 molecule. The shaded areas are used to 
mark orbitals which are localized on the adsorbates. 



the number of adsorbed molecules . The HOMO-LUMO 
gap is almost unchanged except for the cases when the 
orbitals are localized on the adsorbed molecules as indi- 
cated in Fig. [4] by the shaded areas. The orbital energies 
are shifted when the number of molecules increases in- 
dicating the p-doping of graphene. Degradation of the 
bands is clearly consistent with the alteration of the 
charge transfer from graphene to the adsorbates (see Qt r 
in Table QJ. Thus for the adsorption of water, the shift 
of the orbital energy gradually increases when the num- 
ber of adsorbed molecules grows, thereby providing an 
opportunity to control the doping of graphene through 
adsorption of the polar molecules. 



IV. CONCLUSION 

We have investigated the alteration of the electronic 
properties and spin distribution induced by adsorption of 
the polar molecules in nanoscale graphene with oxidized 
zigzag edges. The polar molecules interact with oxygen 
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atoms belonging to the polar C-0 bonds at the edges of 
graphene. When the adsorbed molecule contains hydro- 
gen atoms they interact with the C-0 bonds of graphene 
through generation of the hydrogen bonds (H • • • O-C). 
Several adsorbed H2O or NH3 molecules tends to create 
a cluster due to the strong intermolecular interaction be- 
tween them, but the interaction with graphene aligns the 
cluster structure along the zigzag edge. Graphene with 
oxidized edges tends to donate a charge to the adsor- 
bates and its efficiency is defined by the intermolecular 
distance and by the location of the adsorbed molecules 
relative to the plane of graphene. In the case of H2O 
adsorption, the charge transfer efficiency and the corre- 
sponding shift of the HOMO and LUMO bands as a result 
of p-doping of graphene increase gradually with the num- 
ber of molecules. For the NH3 molecules the maximum 



efficiency of doping occurs when the number of adsorbed 
molecules along the zigzag edge is limited by two or three 
as a result of their closer location to the plane of graphene 
than when more molecules are adsorbed. Therefore, we 
can conclude that when the graphene surface or the edges 
are oxidized, the adsorption of the polar molecules is the 
most viable option for p-doping of graphene. 
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